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Abstract 

A quantum particle moving in a uniform magnetic field and periodic potential (the 
Hofstadter model) has an energy band structure which varies in a discontinuous fashion as 
a function of the magnetic flux per lattice unit cell. In a real system, randomness of various 
kinds should "smooth out" this behavior in some way. To explore how this happens, we 
have studied the dissipative quantum mechanics of the Hofstadter model. We find, by 
virtue of a duality in a two-dimensional space parametrized by the dissipation constant 
and the magnetic field strength, that there are an infinite number of phase transition lines, 
whose density grows without limit as the dissipation goes to zero and the model reduces to 
the original Hofstadter model. The measurable quantity of greatest interest, the mobility, 
can be determined exactly in most of parameter space. The critical theory on the phase 
transition lines has yet to be characterized in any detail, but it has reparametrization 
invariance and defines a set of nontrivial backgrounds for open string theory. 
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1. Hofstadter Model and Dissipation: Introduction and Issues 

The quantum mechanics of an electron moving in two dimensions subject to a uniform 
magnetic field and a periodic potential (which we will call the Hofstadter problem for 
short) has a remarkable fractal energy band structure [Q which varies discontinuously as a 
function of the number of flux quanta per lattice unit cell. Such behavior must, however, 
be "smoothed out" by the inevitable randomness of a real physical system. An interesting 
model of dynamic randomness (due ultimately to weak coupling of the electron to lattice 
vibrations) is provided by the dissipative quantum mechanics of Caldeira and Leggett 
(DQM) P) and the topic of this paper is the application of DQM to the Hofstadter problem. 
A key point is that DQM converts quantum mechanics into one- dimensional statistical 
mechanics (the one dimension being Euclidean time on the electron's world line) with 
long-range interactions strong enough to support phase transitions between, as it turns 
out, localized and delocalized long-time behavior. Thus, the discontinuities of the original 
quantum mechanics are replaced by phase transitions and the problem is to map out the 
phase diagram (in the magnetic field-dissipation constant plane). Presumably the density 
of phase transitions grows as dissipation decreases, recreating discontinuous dependence 
on the magnetic field in the zero-dissipation limit. We are able to show, using a duality 
argument, that this is exactly what happens and that the phase diagram is what was 
found years ago |^ in the study of theta-dependence of the confinement-deconfinement 
transition in Abelian lattice gauge theories! Away from the critical lines, the behavior 
of the system for most of the parameter space is very simple: the long-time mobility, for 
example, is given by the trivial theory for some value of the magnetic field and no potential. 
The critical theories, however, are definitely non-trivial and we have only been able to 
extract limited information about them. Apart from their condensed matter interest, they 
represent new solutions of open string theory . Our results at least indicate that there is 
a rich vein of one-dimensional critical theories to explore, once the appropriate tools have 
been developed. 

The outline of this paper is as follows. Section 2 contains a review of the dissipative 
Hofstadter model. In Section 3 we sketch the arguments showing that, at specific critical 
points, the theory is at a renormalization group fixed point and the mobility can be calcu- 
lated exactly. In the fourth section, we derive the approximate duality transformations of 
the dissipative Hofstadter model, and, in the following section, we demonstrate that the 
duality transformations are exact for a discrete version of the model. In section 6, we use 
these transformations to map out the phase diagram of the dissipative Hofstadter model. 
In the final section, we present a brief discussion of our results. 

2. Setup of the Dissipative Hofstadter Model 

We begin with a brief outline of dissipative quantum mechanics. For details the reader 
is referred to and . The quantum mechanics of a particle subject to a scalar potential 
V and vector potential A is described by the Euclidean path integral 




(2.1) 



with 
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If the quantum variable, X*, is in some sense macroscopic (a good example is the trapped 
magnetic flux in a Josephson junction) there will typically exist an infinite set of degrees 
of freedom to which X'^ is at least weakly coupled and which give rise to dissipative effects 
on the motion of X*. In the classical limit, dissipation can be described by adding a 
phenomenological friction term —rjX^ to the equation of motion. It is then natural to ask 
whether there is a correspondingly simple and universal way to express the quantum effects 
of weak dissipation. A practical motivation is to assess whether coherent quantum effects, 
such as tunneling, can really be observed in macroscopic systems. 

Caldeira and Leggett addressed this problem via a simple model: In addition to 
the coordinates X*, let there be a bath of harmonic oscillators, qa, with a distribution 
of frequencies, Wq,, coupled linearly to the X* with weak coupling strengths, Cq,. If the 
parameters satisfy the functional condition 



E 



^ 4(^-^„) = ^ , (2.3) 



then, when the oscillators are 'integrated out' of the classical equations of motion for X% 
they supply the canonical — ryX* friction term. Since the dependence of the action on the q^ 
is only quadratic, they can be explicitly integrated out of the quantum path integral also. 
The result is a new path integral over the X'^ alone, where the quantum effect of friction 
is contained in a non-local term whose strength is set by the classical friction constant rj: 

ZJ^ = j [VX{t)] e-SQM[x]-s^[x] ^2.4) 

where 



Because of the non-locality of the ry-term, the path integral ( p.4| ) is effectively that of 
a one-dimensional statistical system with long-range interactions. Such systems, unlike 
one-dimensional local systems, have phase transitions (the classic example being the Ising 
chain with interactions ). In the DQM context, the phase transitions are between 
different regimes of long-time behavior of Green's functions (typically between localized 
and delocalized). The crucial qualitative information about DQM therefore concerns the 
phase structure of its long-time behavior and we will analyze the dissipative Hofstadter 
model from that point of view. 

To study the Hofstadter model, we specialize ( p.2|) as follows: X = (x, y) is taken to 
be two-dimensional, the potential is taken to be periodic with period a in both directions, 
and of strength Vq, so that 

V{x, y) = Vo cos{2iTx/a) + Vq cos{2ny/a) ; (2.6) 

and the magnetic field is taken to be uniform and in a linear gauge: 

{A^,Ay) = i{By,-Bx) . (2.7) 
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Leaving aside the potential term, the action ( pT2|) is Gaussian: 



The generating functional for for the full theory is then given by 

Z[r],B,V,F] = J [DX{t)]e~^'^^^~^ (2.9) 

where ^SpfX] is a linear source term, 

Sf[X] = J F{t) ■ X{t)dt. (2.10) 

In order to determine the phase structure of the dissipative Hofstadter model, we will be 

concentrating on the properties of the partition function and two-point function, which 
are obtained from Z[r], B,V, F] as follows: 

Z[rj,B,V]=Z[rj,B,V,0] (2.11) 

and 



(2.12) 



F=0 



Z[77,S,y,0] 5F-{ti)6F^^it2] 

It is useful to define the dimensionless parameters a = r]a /2Tin and /3 = eBa'^/27inc, 
to rescale X by 1/a, and to reexpress the Gaussian action in Fourier space: 

So = ^ f ^{{-^H + ^^")S^J + '^^P^^M X:{u;)X,{u;) . (2.13) 
J ATT Ztv In 

One sees that the dissipation and magnetic field terms are dimension one operators (they 
grow like the first power of uj) and therefore marginal from the renormalization group 
point of view. On the other hand, the ordinary kinetic energy term is dimension two 
and therefore irrelevant: it just acts as a regulator for short times and has no effect on 
universal critical properties. We will set it to zero from now on and use a more convenient 
short-distance cutoff wherever needed. 

The action (127131) with M = plays a very important role: it is the Gaussian fixed 



point theory which governs the long-time behavior of the system whenever the potential 
term ( |2.6| ) is irrelevant. We will shortly identify the region in the a-(3 plane where this 
condition is met. Let us first summarize the essentials of the Gaussian theory: all infor- 
mation is contained in the coordinate two-point function {Xi{t)Xj{t')), the inverse of the 
quadratic form in ( |2.13|) . We will call it Dij{t — t'; z), using the complex number z = a + i(3 
to summarize the parameters on which it depends. Explicitly, 

D,j{t - t'- z) = --^^ log(t - t'f - i-0^sign{t - t')e,j . (2.14) 
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The logarithmic growth of D{t) clearly indicates that the Gaussian system is on the bor- 
derline between localization and delocalization. The coefficient of the 5ij term is called, 
for obvious reasons, the mobility (/i). For zero magnetic field, oo as dissipation goes 
to zero, while for finite field, the mobility goes to zero in the same limit since the parti- 
cles become stuck in Landau levels. The coefficient of the eij term is essentially a Hall 
coefficient: it measures the response transverse to an applied field. 
In Fourier space, the two-point function is 

D,j{u;z) = j^M,j{uj;-), (2.15) 

\UJ\ z 

where 

Mij{uj,x + iy) = xSij + y—eij. (2.16) 

to 

In Section 6, we will make use of the following simple rules for the addition and multipli- 
cation of Mij-. 

M,j{uj; zi) + M,j{uj] zi) = M,j{uj; Zi + Z2); (2.17) 
Mij{uj;zi)M,j{uj]Z2) = M,j{uj- Z1Z2); (2.18) 

and 

cM,j{u-zi) = M,j{u-czi) (2.19) 

for all zi,Z2 & C and c & R. 

Given the explicit Gaussian two-point function, one can easily calculate the effective 
dimensionality of a perturbation such as ( |2.6D and identify the regions in the z-plane where 
it is irrelevant. This is a simple one- loop renormalization group calculation and it is easy 
to show that (|2.6| ) is marginal when the mobility coefficient a/(a^-f/3^) is unity and 
irrelevant when a/{a'^ + /3^) is greater than one. This condition defines a circle in the 
2;-plane, of radius one-half and centered at z — (1,0), inside of which (|2.6D is irrelevant 



and on which it is marginal. For zero field, the critical value is a = 1, a fact which has 
been known for a long time . Outside the circle, ( p.6|) is relevant and, if one knew that 
there were only two phases, one would have to identify that entire region with a localized 
phase. This is almost certainly correct on the zero magnetic field line, but, given the rich 
dual and self-similar structure of the pure quantum mechanics of the Hofstadter problem, 
certainly wrong for non-zero field. 

Two loosely related problems arise at this point. The first is to map out the phase 
structure in the region outside the "first circle", where the periodic potential is relevant. 
We believe we have a complete solution of this problem. As will be explained in subsequent 
sections, the key to the solution is an approximate duality symmetry (generated by two 
transformations z —* 1/z and z z + i) which relates behavior in different regions of the 
2-plane. The second problem is to characterize the family of critical theories which live 
on the boundary of the first circle (and its images under duality). At present we know 
very little about them, although, as we shall explain in the next section, they become 
particularly simple at the points where P/a is integer. There is a further motivation for 
studying these critical theories, whose elaboration must await another paper. It is that, 
as explained in ^ , these critical theories satisfy an infinite set of Ward identities deriving 
from reparametrization invariance, possess a manifest SL{2, R) symmetry, and represent 
nontrivial solutions of open string theory! 
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3. Fermionization and Exact Fixed Points 



In this section, we summarize some of the properties of the critical theories at the 
points on the "first circle" with integral Using a convenient regulator, we show that at 
these points the theory has no logarithmic divergences and is therefore at a renormalization 
group fixed point. In addition, we show that the theory can be fermionized and that the 
free energy and mobility can be exactly calculated. 

In order to calculate the partition function and correlation functions of the X^^s, it 
is most convenient to expand the path integral in powers of the cosine potential. We 
then have to evaluate correlators of products of X^'s and e^'^^ 's (the latter coming from 
expanding the cosines as a sum of exponentials) and the result can be written in terms of 
integrals over products of factors of the type {X^^(ti)X^{tj))Q and e^'^-^^^^^^^^^^^^^'o^ where 
the expectation value is with respect to the gaussian propagator Dij{t — t';z) and the 
integrals are over the t coordinates of the potential insertions. In order to get finite results 
at all stages of the calculation, we need both an infrared and an ultraviolet cutoff. We 
obtain the former by putting the parameter t on a circle of circumference T, and the latter 
by multiplying the Fourier-transformed propagator Dij{uj] z) by e"'^''^' (we do this instead 



of giving the particle an explicit mass). With these regulators ( |2.14| ) is no longer correct, 
but it is not too hard to show that 



(X^(ti)X^(t2))o = -h7r 



27rz 



a 



+ 



w\ — W2 



T J a"^ + P"^ \_{wi — e ^W2){wi 



1 + 



)WiW2 



{w-i 



'W2){wi 



e'=W2 



(3.1) 



e W1W2 



{wi - e''w2){wi - e ^W2) 



X 



wi — W2e 
W2 — wie'^ 



W2 



1 ±e 



(3.2) 



where we have defined a new position coordinate Wj = e^'^**^/^ which lies on the unit 
circle. Under this transformation, integrals over tj are converted into contour integrals 
of Wn around the unit circle according to f dt ^ S We note that whenever 

ajip?' + /3^) = 1 and /3/ct is an integer, the expressions in equations (|3.1[) and (|3.2|) are 



rational functions of Wj and e*^, so the integrands for the partition function and n-point 
functions will also be rational functions. Furthermore, because the integrand factors into 
a product of simple poles, we can perform the integrals by successively extracting the 
residues of the poles in the integration variables which lie within the unit circle. At each 
stage we are left with a rational function of the remaining Wj's and e*^. This implies 
that, in the end, we can only get pole divergences in e as e — > 0, and no logarithmic 

divergences. In particular, for correlators of marginal operators such as X{t\ which can at 
worst have logarithmic divergences, there can be no divergences at all. As a result, when 
a I (a^ -|-/9^) = 1 and [3/a E the theory should be at a zero of the /9-function. As we will 
show elsewhere, this also implies that the theory satisfies an infinite set of Ward Identities, 
which derive from reparametrization invariance. (As was explained in 0, these Ward 
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identities are what make possible a string theory interpretation of systems of this kind.) 
This is very much hke the 2D case, where scale invariance implies conformal invariance. 
We note that the absence of logarithmic divergences occurs for any value of the potential, 
Vq, which means that these fixed-point theories are really fixed lines. 

We can use the regulator and contour integration procedure just described to calculate 
the partition function and n-point functions to any order in Vq, but it becomes extremely 
tedious for higher orders. At the special point a = 1, = 0, we can obtain a number of 
results exact to all orders in Vq by fermionizing the unregulated theory and then regulating 
the fermionic theory. We observe that if we express the X's in terms of two free fermions 
"tp+{t) and i{^-{t) according to 

e^^W=V.;(t)^_(t), (3.3) 



and 



2ni 



X{t) = — ^lm+{t) - i^im-{t) . (3.4) 



and define the fermion propagators to be 

{i^l{t^)Mt2)) = (^l(ti)^_(t2)) = '-^^^ (3.5) 

W2 — Wi 

and 

(V^i(ti)^_(t2)) = (^l(ti)^+(t2)) = , (3.6) 

we reproduce all the bosonic correlators of X's and e^*^. Note that the set of operators 
under discussion here generates an SU (2) algebra and is a particular example of the discrete 
states which have been extensively discussed recently in two-dimensional quantum gravity 

Since the potential term, and all the operators of interest, are now quadratic in free 
fermions, the unregulated problem reduces to free field theory. The ultraviolet-regulated 
theory might not, however, and care must be taken here. We regulate by multiplying the 
Fourier transform of the fermi propagators by e"'^''^', and again find that all the integrals for 
the free energy and n-point functions are contour integrals of factored rational polynomials. 
This time, all the integrands are just products of loops of the form ^1^=1 l/('^i ~ e~'^t(;i_|_i), 
with wi-^-i = wi, so we can easily perform the integrals at any order in Vq. Furthermore, at 
least for the partition function and the two-point function, it is easy to resum the expansion 
in powers of the potential. The result for the free energy density is 



11 /T/ T \ 

Y ^ i^^-l)'' (M:e-(^(0)x(0))^ +o(e), (3.7) 



N=l 

while the result for the two point function is 

Note that both of these expressions depend on the peculiar dimensionless parameter 
VQTe~^'^^^^'^^^^^ . In the limit that the ratio of infrared to ultraviolet cutoff is large. 
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this parameter reduces to Vo/e = Vq. The only renormahzation needed to render the two- 
point function (and, indeed, aU n-point functions of the X's) finite is to choose Vq to be 
finite. This doesn't quite render the free energy density finite: it has a 1/e piece which we 
must subtract away by hand. We only attach physical significance to the finite piece left 
over after this subtraction and we note that in this particular case [a = 1 and /3 = 0) this 
finite piece is zero. Note also that the functional form of the two-point function follows 
from SU{1, 1) invar iance if X transforms with weight unity. Both of these points and their 
string theory significance are explained at greater length in [|5[. 

Apart from an overall Vo-dependent, finite, multiplicative renormahzation, ( p.8|) is 
the same as the Gaussian two point function D{t — t']l) transformed to the unit circle 
(consider (|3.1| ) at a = 1, /3 = 0). So, at the a — 1 fixed point, the ratio of the renormalized 
mobility to its Gaussian value has a specific dependence on the potential strength: 

fieff/fiGauss = (1 + (e-<^(°)^(o)VoT/2)2)-i = (1 + (^0/2)')-^ . (3.9) 

This is in rough accord with the results of Guinea et. al. Since they worked in the 
tight-binding approximation, rather than in the weak-potential expansion used here, the 
renormahzation of the Gaussian two-point function is different in the two cases. Rather 
surprisingly, even for non-vanishing potential strength, the higher n-point functions turn 
out to vanish for non-coincident points (i.e. when ti 7^ tj for 1 < i, J < n). On the other 
hand, when all the ti are equal, the regulated 2n-point functions are proportional to 
Thus, the n-point functions for n > 2 are zero except for contact terms. Because of the 
rescaling of the 2-point function and the contact terms for higher n-point functions, the 
fixed point theory when a = 1 does not appear to be a simple gaussian theory. Precisely 
what it is is not yet known in any detail. 

It is also possible to fermionize the theory at the other special points where G 
Z and a = + Although the X's require more care to fermionize because 

(e*^*^*i^X(t2)) oc ^sign(ti — ^2) is ill-defined without the regulator, the results for the 
free energy and 2-point functions are very similar to what we found at the a = 1 fixed 
point. 



4. Duality Symmetries from Coulomb Gas and Instanton Methods 

We will now show that the dissipative Hofstadter model has at least an approximate 
duality symmetry under the modular group consisting of the transformations z —>■ 1/z, 
z ^ z + i and all their compositions, where z = a + i(3. In Section 6 we will use this 
duality to obtain an overall picture of the phase structure of the theory. To demonstrate 
the symmetry, we begin by deriving a generalized Coulomb gas representation of the DQM 
path integral. The first duality symmetry then follows from a manifest invariance of the 
Coulomb gas. Next, we derive a second Coulomb gas representation, which is related 
to the first by a transformation of z and Vq. Insofar as both Coulomb gases are valid 
representations of the same model, we obtain the second duality transformation. These 
two transformations on z are enough to generate the full modular group. This strategy was 



first applied to DQM without a magnetic field by Schmid [|T0| to obtain restricted duality 



under the transformation o; — > 1/a of the dissipation constant only. It was also used in 



reference [11 for the Hofstadter model with no dissipation. 
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For the first derivation of tfie Coulomb gas, we simply expand the path integral in 
powers of the potential and further expand the cosine potential as a sum of positive and 
negative frequency exponentials: cosX(rj) = h'^e-=±i^^^^^^'^'^ ■ Because the potential 
( pj.6| ) is a sum of terms periodic in X and Y, we must actually define a two-dimensional 
vector of charges 6^. A collection of n of these charges, located at times r^, can be associated 
with a charge density 



3 = 1 

With these definitions, the generating function (|2.9|) becomes 

Z[z, Vo, F] = l [DX{t)] f^jdn... dr^ (y ) " ^ E 



(4.1) 



e ^"i. 



(±1,0) 
J~ (0,±1) 



where 



Sq = Sg + SF-i / X{t) ■ p^{T)dr. 



(4.2) 



(4.3) 



Because Sq is quadratic, the X integration can be done exactly, with the result 



Z[z,Vo,F] =CZ, 



gen 



Sgen{z, F); 



Vo 



where the action, Sgen{z, F), is given by 
1 



/I \ 
dri... dT2Nj^^^ y-jj exp[->Sgen(2, F)] 



(4.4) 



Sgen{z,F) = ^ / drdr' 



iF{T)+f{T) ' D{t-t'-z) iF{T')+p^{T" 



(4.5) 



The allowed values for the vector "charges" in the expression for Zg^^ are (±1,0) and 
(0, ±1). Integrating over the zero mode of X gives the additional requirement of charge 
neutrality: ^j^^ = 0. 

To obtain the partition function, we set F = in equation (4^). It will be useful 
to note that the resulting object is identical to the grand canonical partition function for 
a certain generalized Coulomb gas. Let us define a Coulomb gas interaction "energy" 
between a collection of charges at locations Ti by 



Scg{A + iB) = l^e,-I)(r,-r,;— i 



A + iB' 



-■3 1 



(4.6) 



where D is the Gaussian two-point function of ( 2.14 ). Writing out the full expression for 
D, we obtain the explicit expression 

Scg{A + iB)=AY,e,- e, \n{U - t,f + J] sign(r, - r.^e'^^ef e^". (4.7) 



i<3 
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This is just a collection of two-body interaction terms, with real and imaginary parts whose 
strengths are set by A and B respectively. The real part looks like a standard Coulomb 
potential term. The imaginary part is of course unconventional and that is why we refer 
to this as a generalized Coulomb gas. Its partition function is 

ZcG[A + iBX]=J2Yl I (^ri...(ir2Ar-l-y(C)'^exp[-^CG(A + zi?)], (4.8) 

where C is the fugacity and the sum is over even numbers of charges only because we impose 
overall charge neutrality. The imaginary part of the interaction energy gives a phase to 
each term in the partition sum. In this regard it is similar to the 6 term of QCD. 

The partition function of interest to us, (|4.4|) evaluated at F = 0, can thus be written 



in terms of a generalized Coulomb gas as follows: 

Z[z,Vo\=CZcg[1/z,Vo/2\. (4.9) 

In short, the original system is equivalent to a generalized Coulomb gas in which each 
particle has charge Cj = (±1, 0) or (0, ±1) and fugacity C, = Vo/2, and the bonds between 
particles are assigned both an energy and a phase determined by the parameter z ~ 
l/{A + iB). 

By varying the generating function ( [4.4| ) with respect to F and setting F = one 
gets expressions for n-point functions of the X's in terms of Coulomb gas correlators of 
the charge density p. For the two-point function, in particular, we obtain 

{X^{uj)X-{-uo)){z, Fo) = D^'^iu:- z) - D^'^iuj; z) {p%u:)~p\-u:)){l/ z, Vo/2) D>^^{uo- z), 

(4.10) 

where by (p'^(u;)p^(— Ci;))(l/2;, V()/2) we mean the density-density correlation function in 
the Coulomb gas with coupling parameter Xjz and fugacity Vb/2. We should emphasize 
that although we have expanded in powers of Vq, no terms have been dropped, so our 
results to this point are exact. 

The generalized Coulomb gas described by ZcG[A+iB, Q possesses a symmetry under 
shifts in A+iB. Because all the charges are integral, the phase factor appearing in a generic 
term in the partition sum, 

exp[z7ri? J2 sign(^. " ^jOe'^'ef e,"] , (4.11) 

is invariant under B B + 1. Consequently, both the partition function of the generalized 
Coulomb gas and the charge-density correlation functions remain invariant under the shifts 
A + iB ^ {A + iB) -f m, for integer n, so that 

Zcg[A + iB + in, C] = Zcg[A + iB, Q, (4.12) 

and 

{p''{u)p\-u)){A + iB + inX) = {p''{uj)p\-u)){A + iBX). (4.13) 
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If we let 1/z = A + iB and 1/z = A + iB + in, then we can use these equations and 
equations ( |4.9| ) and ( ^.101 ) to solve for the transformation of the partition function and 
two-point function of the Hofstadter model when z goes to z — z/{l + inz). We find that 
the partition function remains the same and the coordinate two-point function transforms 



as 



{Xiu)Xi-u)){z,Vo) 



D{uj; z )+uj^D{uj; 1 + inz) {X{lv)X{-uj)){ 



n 



1 + inz 



Vq) D{uj; 1 + inz) 



(4.14) 

We believe this transformation to be exact because we have summed to all orders in 
perturbation theory. There are some subtleties having to do with regulation, but they 
should have no infiuence on the long-time behavior. In the next section, we will describe 
a discretized version of the dissipative Hofstadter model for which duality transformations 
like ( [1.14| ) are exact without any qualifications. 

The second Coulomb gas is obtained by making an instanton expansion of the path 
integral. Here we assume that if the friction is small enough, the main contribution to the 
path integral comes from the "approximate" solutions to the Euclidean classical equations 
of motion in the absence of friction. In these solutions the particle moves between adjacent 
mimina of the potential along instanton paths, which correspond to the particle tunneling 
from one minimum to another. These frictionless instantons for the Hofstadter model were 



studied in some detail in and we will draw heavily on results presented there. When 

M = 0, the classical action for the instanton is given by s = e ^^'^ with 7 ^ 7.33, and 
when /3 ~ 0, it is s = 4\^VoMa^/nh. When M = 0, the instanton solution starting at 



(x, y) = (0, 0) and ending at 2nej, with ej = (±1, 0) or (0, ±1), is given by | |TI[ ] 



where 



and 



/(r) - ejh{T) + {zx ej)g{T) 
h{T) = 2tan"^ 



2V0 



sinh(47r~ — -r) 
ah 



(4.15) 
(4.16) 



g{T) = —i cosh ^ 



1 -f2sech(87r2^r) 
ah 



(4.17) 



We can add together a series of these "kinks" to construct the approximate solution 



^n{r) = ejh{T -Tj) + {zx ej)g{r - Tj), 



(4.18) 



whose Fourier transform can be written in terms of the charge distribution defined in 
equation ( [1.1| ): 

X^(w) = h^p^ +gcuizx p^). (4.19) 

When we substitute this solution back into the full action and then sum over all such paths, 
we once again obtain a Coulomb gas expression for the partition function: the two signs 
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of the charge correspond to instantons and anti-instantons; the locations of the charges 
correspond to the instanton center of mass collective coordinate; the fugacity is given by 
the exponential of minus the instanton action times a fluctuation determinant {Ke~^), 
where K = Vo/>/\P\ when a = and /?>>!; and the instantons acquire long-range 
Coulomb interactions and phases from the dissipation and magnetic field terms in the 
action. The expression for the generating functional has the same structure as before. It 
is 

Z[z,Vo,F] = Zg,^[S,^,t{z,Fy,Ke-'], (4.20) 



where now the energy function, .Sgen, in equation (|4.4|) is replaced by 

jk 

Here, A specifies the interaction energy and phases between two instantons and it is similar 
in structure to the Gaussian interaction function D{Tj — Tk;l/z). In particular, its off- 
diagonal terms are exactly the same, but the diagonal terms depend in detail on the shape 
of the instanton solution. However, in the long-time limit \{Tj — Tk)4:iT'^Vo/{ah)\ >> 1, A 
turns out to have a universal form 

A{tj - Tk-, z) ~ D{tj - Tk] l/z). (4.22) 

Thus, in this limit the instanton interactions are identical to those of a Coulomb gas with 
interaction parameter z instead of l/z. We will assume that properties of the critical points 
of the theory are sensitive only to this long-time limit. In that case, the partition function 
can again be written in terms of the Coulomb gas partition function given by equations 
(U) and (13). It is 

Z{z, Vo) = NZcg{z, K{Vo)e-'), (4.23) 

where denotes the value of the partition function in the absence of tunneling. 

We can again derive an expression for the two-point function in terms of the charge 
correlation function of the new Coulomb gas: 

{X^{u)X-{-u)){z, Vo) = H^^u) {~p%u)~p\-u)){z, Ke-n H^^iu;), (4.24) 

where -ff(w) is 

H^"'{uJ) = h{uj)S'"' + g{u)e^r (4.25) 
The small u (large r) behavior of h and g is such that ( ^.24[ ) reduces, as u; — 0, to 



{X^{u)X'^{-u)){z, Vo) = D^'^iu; 1) {p^{u;)p\-u)){z, Ke"^) D^''{u;; 1). (4.26) 

The instanton method should be valid whenever the classical action of the Hofstadter 
model (without dissipation) is large and the coefficient of friction is much smaller than 
the classical action. We also assume that the the large-time behavior of the system is not 
affected by the particular (regulator-dependent) small-distance form of the propagators. 

To summarize, the preceding calculations demonstrate that the original partition func- 
tion is equivalent to the partition function for two different Coulomb gases of "X" and "y" 
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particles which correspond to the X and Y components of e'j, respectively. The particles 
have charges ±1 and a fugacity which depends on which Coulomb gas is under consider- 
ation. The structure of the interaction energy is the same for both gases: like particles 
(both X or both Y) interact via — 2Aefe^ log \ Ti — Tj \ for \Ti — Tj\ » 1, and unlike particles 
via Z7ri?sign(ri — Tj)e'^'^e^ej . When we treat the potential as a perturbation, the fugacity 
is given by Vo/2 and the interaction strengths hy A + iB = 1/ z (where z is the parameter 
introduced in ( |2.14|) ). When we use the instanton method, the system is described as a 
gas with fugacity / zeta — K{Vo, z)e~^ and interaction strength A + iB = z. Therefore, 
we have established that the system with dissipation and flux per unit cell given by z and 
with a potential of strength Vq is dual to the system with z ^ 1/z and a new potential 
Vo/2 = K{Vo,z)e-'. 

The two-point function of the original system can be expressed in terms of the charge- 
density correlation function for either of the two Coulomb gases described above. Upon 
eliminating the charge density correlators between the two expressions, we get the follow- 
ing duality relation between two-point functions of the original theory in dual regions of 
parameter space (we expect it to be valid for large times, or, equivalently, small iv): 

{X{uj)X{-uj)){z, Vo) = D{uj] z) - u?I){uj; z) {X{uj)X{-uj)){\I z, V^) D{uj; z) . (4.27) 

Note that the fugacity transforms along with z [Vq —>■ Vq). This sort of relation has 
the strongest consequences in circumstances where there is in fact no dependence of the 
two-point function on fugacity. As we shall see, this covers a lot of territory. 

Lastly, we may combine the "exact" duality transformation z — > z/{l + inz) with 
the approximate duality transformation z ^ 1/z. A complication here is that under 
the second duality, Vq transforms in a 2;-dependent way. Consequently we learn that the 
dissipative Hofstadter system at z and Vq is dual to the system ed, z + ni and some other 
potential strength. The two systems should have the same partition function up to an 
analytic prefactor. Additionally, if for some reason the two-point function is not explicitly 
dependent on V, then, by composing ( [4.27| ) and ( ^4.14|) , we obtain 

{X{uj)X{-uj)){z + in) = {X{uj)X{-uj)){z) . (4.28) 

We conclude that the system has an exact symmetry under z ^ z/{l + inz) and an 
approximate symmetry under all compositions of z ^ 1/z and z — * z + i. We will draw 
the consequences of this for the phase structure of the theory in the final section. 

5. The Discrete Gaussian Model: an Instructive Example 

Because our model has an approximate duality under z ^ z + i and z 1/z, it is 
natural to suspect that there is a similar model which has an exact duality under these 
transformations. The discrete gaussian chain with the interaction V{r) = otjir'^ + 1/4) is 
known to be exactly dual to one with a replaced by (1/a) . We find that if we couple 
two such chains via a magnetic interaction, the system is then exactly dual under z ^ z^i 
and z -^Xj z. We define such a system as follows. There is a one-dimensional lattice located 
at the sites z = 1, 2, and a second dual lattice at i' = 1/2, 3/2, A— 1/2. At each site 
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there are height variables and h^, with /if_^jv — h^^^N ~ ^1'- '^^^ height variables 

can take on any integer value except for and h^j^_i/2^ which we set to zero in order to 
keep the partition function finite. We will define our discrete Gaussian Hamiltonian with 
magnetic interaction to be 



-l2n^f3Y.hUhl,-hy_,) + l2n^pJ2hm'+^-h^'-^: 



(5.1) 



where V is chosen so that its Fourier transform V satisfies 

k 

WN{k) = Vn{0) - VN{k) = 27r|sin-|. (5.2) 

The partition function 

^^^w = EE^"''^^ (5-3) 

hf h^, 

depends on two parameters a and 13, which we have combined into a complex parameter 
z = a i(3. With these definitions, the momentum space version of ( |5.1[ ) reads 

H^"^ = ^J2^.k-D-\k;z)-hk (5.4) 

k 

where and are the Fourier transforms of /i| and , 

_ ^e^^^Ti^, hi = Y^e'^^'hl, (5.5) 

k k 

(where k = 2'Kr/N and r runs over the integers satisfying —N/2 < j < N/2), and 
D-^{k;z) = 2WN{k)M{k;z), with M defined by 

a 3— 

M{k-z)=\ I'^l 1. (5.6) 



a 



D{k; z) is the discrete-time version of the propagator D{u);z) originally introduced in 

Now we can do a Poisson resummation trick, exactly as in [^, to reexpress the 
discrete Gaussian system as a Coulomb gas. The procedure is to replace the hi and hi/ by 
continuous variables Ui and Ui' (over which we will integrate) but to include delta functions 
of the form 6^^*'^^'^^ to pick out the integer values of u. On integrating out the 

Ui^s we obtain 

Z^%] = NC^{a^ + /32)-(^-i)/2Z^«[l/^], (5.7) 

13 



where = Ylk^o^/^i^) is given by 



„x y 



7r2 



(5.8) 



Here, qk = {q^-, q^), where q^ and q^ are the Fourier transforms of qf and q^,, respectively. 
The qi are constrained to satisfy — — ^- not too hard to recognize ( |57^ ) 

as a discretized generahzed Coulomb gas. Indeed, in the limit oo, the interactions 

between the charges are exactly the same as in the model of the previous section. The key 
differences are that: the charges can now take on any integer value (they were previously 
restricted to ±1); there is only one charge per lattice site; and the fugacity is 1 rather than 
being freely variable. 

Next, we perform Cardy's transformation on the expression for Z^'^ . We define 
the new height variables. If, and , by 

Qi = ^i+l/2 ~ Qi' = ^i' + l/2 ~ ^^'-1/2' ^0 = ^1/2 = 0- (5-9) 

On making this variable substitution in ( |5.8|) we find, on comparing with ( |5.4| ) , that we 
have reproduced Z^*-^ with z 1/z. Specifically, 

Z^^[z] = NCMia' + P')-^''-'^/'Z^^[l/z]. (5.10) 

This exact duality between the system at z and 1/z depends on the precise form (|5.2|) and 
( p.6| ) of the interaction. The partition function in equation ( |5.3| ) is also manifestly invariant 
under z z + i. Thus, the system is exactly dual to itself under all compositions of the 
two transformations z —>■ 1/z and z —>■ z + i. We will shortly explore the consequences of 
this. 



By a slight variant of the above arguments, one can derive, following [|T2| , an exact 
duality relation for the two-point function: 

{h{kyh{-kY){z) = D^^'ik; z) - M^P{k] 1/z) (h{k)Ph{-kf){l/z) M^^{k- 1/z). (5.11) 

The duality relation ( |4.27|) for the continuum problem is equivalent to this if, for some 
reason, the two-point functions don't depend on the potential amplitude or there is some 
self-dual value of the fugacity for which V = V. 

At the self-dual point, 2 = 1, equation ( |5.11J ) determines the two point function 
exactly. It is {h{k)^h(—k)^){l) = ^i)(/c; 1), which is half the value for the continuous 
gaussian model. (Note that there is a certain discontinuity here: anywhere inside the first 
circle, no matter how close to z = 1, the two-point function is given by D{k, z).) We can 
apply the transformations z ^ z + i and z —>■ 1/z to this expression to find the value of 
the mobility for all z given by 

d — ib 1 .ab + cd ^ , , ^ , , ^ 

z = = — ;5 5" — z— 5 5", for ad — bc = l, a,b,c,dEZ. (5-12) 

a + ic a'^ + + 

It is 

{Kh^_^){z) = \{a' + ^)b{kA). (5.13) 

We observe that at these special points (which include the special points on the first circle 
mentioned in the previous section) the off-diagonal propagator (or the Hall coefficient) 
vanishes even though the magnetic field at these points is not in general zero. 
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6. Duality and the Phase Diagram 

Now we want to use these results to map out the phase structure of the dissipative 
Hofstadter model in the (a, (3) plane. Along the zero-field line, there is no doubt that the 
system is localized for a > 1 and delocalized for a < 1. The Coulomb gas treatment of 
Section 4 allows us to expand on this somewhat. There are two 'dual' Coulomb gases: in the 
first, the charges correspond to insertions of the cosine potential and the fugacity is one- half 
of the potential strength Vq; in the second, the charges correspond to instantons describing 
tunnelings between potential minima and the fugacity is essentially the hopping probability 
of the particle. The one-loop renormalization group argument applied to a Coulomb gas 
says that fugacity scales to zero (and the associated charges become irrelevant) when the 
coefficient in front of the logarithmic self-interaction energy is greater than a critical value. 
Applied to the first Coulomb gas, this argument says that the periodic potential insertions 
are irrelevant everywhere within the circle a'^ + f3'^ = a. Applied to the second Coulomb 
gas, it tells us that the fugacity for instantons scales to zero (i.e. the hopping operator is 
irrelevant and the particle is localized) everywhere to the right of a = 1. Thus, it seems 
very plausible that the system is delocalized within the 'first circle' -|- = a and 
localized in the half-plane to the right of a = 1. As partial confirmation of this, recall 
that, in Section 3, we were able to exactly calculate the two point function at those points 
on the 'first circle' where (3/a & Z and indeed found delocalized behavior. 

With the help of duality arguments, we can parlay these two results into a global 
picture of the phase structure. There are two levels of duality which we could imagine 
using: The discrete Gaussian model has an exact invariance under the two-element modular 
group z ^ z + i and z — as described in Section 5. The Hofstadter model, however, 
is exactly invariant only under the subgroup z z/{l + inz) (for n integer), as described 
in Section 4. What's missing for the Hofstadter model are the transformations z z + ni 
or z ^ 1 / z: under those parameter transformations, the strength of the periodic potential 
transforms as well. However, in the regions where the potential is effectively irrelevant 
(not just inside the first circle, as we shall see), the shifts and inversions should still be 
valid symmetries of the asymptotic behavior of the theory. 

As a first step, let us see how far we can get with the exact symmetry of the Hofstadter 
model. Under z ^ z/{l + inz), the half-plane ct > 1 is mapped into the circle of radius 
l/(2n^) and center l/(2n^) — i/n. This is tangent to the first circle at a point where 
f3/a = —n and also to the /9-axis at /3 = —1/n. These circles are displayed in Fig. 1. Now 
we use the exact relation ( [4.14|) to compute the mobility (i.e. the long-time or low-frequency 
limit of the two-point function) in the circles in terms of the mobility in the region a > 1. 
Since this latter is a region of localization, its mobility is zero, and (|4.14[ ) expresses the 



mobility in the circles in terms of the mobility of the Gaussian theory, encoded in ( |2.14| ) . 
The result is 

(X^mX^miz^V) = D^"'(t;z+ -) for |t| » 0, (6.1) 

n 

which is to say that, as far as mobility is concerned, the n-th circle is just a translation 
parallel to the (3 axis of the Gaussian mobility evaluated within an appropriate sub-circle 
of the first circle. The simple renormalization group argument implies that the mobility 
is indeed Gaussian within the first circle, so this is consistent with the z ^ z + i duality. 
These results are completely different from what one would find if the periodic potential 
^) were zero: in that case the mobility would be given exactly by ( |2.14]) everywhere in 
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the 2;-plane. So, although there is no exphcit dependence on Vq, the potential is clearly 
relevant outside the first circle! 

Because the mobility inside the first circle and all the circles shown in Fig. 1 has no 
explicit F-dependence, we should be able to use the duality relation under z ^ z + in 
(equation ([4.28D ) to find the mobility for all points that are translates of these regions. 
As indicated above, this symmetry is probably correct for the asymptotic behavior of the 
Hofstadter model. This claim is supported by the observation that for the two circles 
immediately above and below the first circle, we have just shown that the mobility is 
D{t] z^i), which is exactly what we would expect from the duality transformations given 
by equation ( [4.28| ). When we apply the translation by im to the first circle, we obtain a 



chain of circles tangent to each other and to the lines a = and a — 1. Next, we can 
repeatedly apply z z/{l + inz) and z ^ z + in to these circles, which starts filling in 
the regions between the line a = and the first chain of circles with new circles that are 
tangent to the /9-axis and already existing circles. The result is shown in fig. 2. It is the 
same phase diagram found by Cardy for abelian lattice gauge theories with theta terms 0] . 
In the diagram, between any two circles and the /3-axis, there is another, smaller circle. So, 
the boundaries of these circles, which, as explained below, should identify lines of phase 
transitions, get denser and denser as one gets closer and closer to the a = line (the 
zero dissipation line). On this line, every point with rational P is tangent to a circle, and 
every point with irrational jS is not. This connects nicely to the Hofstadter model (with 
no dissipation) where the behavior is quite different depending on whether /3 is rational or 
irrational. 

According to the equations ( [4.14|) and ( [4.28| ) for the transformation of the two-point 
function, we find that inside all the circles the mobility is Gaussian: in a circle tangent to 
the /3-axis ai z — p/q, with p and q relatively prime, the mobility is given by D{t; z — ip/q). 
The triangular regions between the circles belong to some other phase which this argument 
doesn't help us identify. 

Along the phase transition circles, one expects to find a non-trivial critical theory. 
Each of these circles is a mapping of the critical line at a = 1 (or the critical circle 
a'^ + P'^ = a) by compositions of the transformations z ^ z + i and z 1/z. Under 
these transformations, the partition function is multiplied by an analytic function (which 
depends on C and in equations ( [4.9|) and (|4.23| )) and the strength of the periodic 
potential changes. Because the existence of the phase transition on the first circle is 
independent of Vq (at least for small Vq), neither of these two changes to the partition 
function should affect the fact that it is critical. Consequently, we expect all the circles 
to be critical circles, at least for some range of values of the potential strength. However, 
because we do not yet know exactly how Vq changes under these transformations, we know 
little about the critical theories, except at the points where P/a = n for n an integer, 
which is where the other circles are tangent to the first circle. For these points, we can 
apply the exact duality transformation ( [4.14|) to the two-point function at z = 1 (given 



in equation (p?8|)) and then simplify the resulting expression by using equations ( |2.17|) , 
( P^) , and (Pl^)to obtain 

{X{t)X{0)){z^, Vo) = D{t- z^) 4- (/(Fo) - 1) B{t- zl) (6.2) 

where z^ = 1/(1 + in) and /(Vq) is the factor by which the two-point function at z = 1 
differs from the Gaussian two-point function ( |2.14]) . With the specific regulator used in 
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Sect. 3, we found that 



fiVo) = (1 + (]^e-<^('^)^W))2)-i = (1 + (^0/2)^)-^ , (6.3) 

but, in general, its exact form depends on how the theory is regulated. 

We note that this two-point function depends on the strength of the potential. For 
the discrete Gaussian model, there is no parameter corresponding to Vq, but the two-point 
function is consistent with equation (|6.2|) with /(Vq) = 1/2 or Vq = 2. For this special 



value of /(Vb), the two-point function at Zn = 1/(1 + in) becomes 

{X{t)X{0)){z^, Vo) = ^(1 + n^)D{t; 1). (6.4) 

A point to emphasize here is that the off-diagonal part of the mobility (the Hall coefficient) 
vanishes at these points for this particular value of Vq, even though the magnetic field is 
not zero. These considerations make it quite plausible that the discrete Gaussian model of 
Sect. 5 has the same critical behavior as the continuum dissipative Hofstadter model for a 
special value of the potential strength. 

In sum, the potential strength is an exactly marginal variable. The higher-point func- 
tions for a = 1 reduce to products of two-point functions plus contact terms (which have 
not been explicitly computed) and for the other points on the first circle we have evidence 
that some of the higher n-point functions are non-zero, SU(1,1) invariant functions. There- 
fore, even if we don't know the details, we are guaranteed that the critical theory is not 
identical to the Gaussian theory 



7. Loose Ends and Problems for the Future 

The main result is that we have discovered a phase structure for the Hofstadter model 
(and a discrete version of it) which shows how the fractal energy level structure of the pure 
quantum mechanics is matched by a fractal structure of phase transitions in the dissipative 
theory. There are several issues to pursue. We have identified critical surfaces with critical 
behavior depending on at least two parameters: potential strength and arc position in the 
a-(3 plane. We know little about the detailed properties of these critical theories and it 
would be very interesting to study them in more detail, especially in the discrete Gaussian 
model, where things should be simple. We have seen in Section 3 that at the special points 
on the first circle something remarkable happens. We have hints that at these values of a 
and /3, most of the n-pt functions (for n > 2) reduce to contact terms, which is reminiscent 
of topological field theory, but we do not yet know how to pursue this connection. These 
critical theories have an interpretation as solutions of open string theory with nontrivial 
backgrounds. It will obviously be interesting to explore the spacetime interpretation of 
these states. We also want to understand the nature of the phases in Fig. 2 which fill in 
the interstitial regions between the circles. They will presumably be some sort of Coulomb 
phase which is neither localized nor delocalized. We have yet to determine what this means 
quantitatively for the behavior of the mobility and Hall coefficient. 

Another line of development takes us into N-body physics. We have seen that the 
transport coefficients (especially the Hall coefficient) of a single electron have remarkable 
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structure in the presence of dissipation. We can actually turn the Hall effect off by choosing 
the right potential amplitude at special values of the field. Filling Landau levels with such 
electrons would doubtless lead to new quantized Hall effect phenomenologies. There is a 
speculative connection with nonperturbative open string theory, because inserting many 
boundaries into the string possibly corresponds to having many electrons in the dissipative 
quantum system. 

The final big issue is the possibility of constructing experimental realizations of the 
physics we have uncovered. Work is already under way to observe, separately, the Hof- 
stadter spectrum and the macroscopic tunneling effects for a two level system. It may 
be possible that specially constructed junction arrays could be used to experimentally 
investigate the combined effects of the magnetic field, periodic potential and dissipation. 
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